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Explaining superfluidity. 




Fritz London suggested using Bose Einstein condensate (BEC) 
to explain superfluidity. This are possible for 2 reasons: 



1. Helium-4 is a boson, so it may undergo Bose Einstein 
condensation. 

2. The liquid would then become a wavefunction, so it would 
show no viscosity unless it flows with enough energy. 
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1. Fritz London's explanation of superfluidity in liquid helium-4 
using Bose Einstein condensate (BEC). 



2. Theory and experimental search for BEC. 



3. Laszlo Tisza's 2-fluid model theory. Experiments. 



4. Lev Landau's critical velocity experiments. 



5. Explanation on perfect conductivity. Second sound 
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In order to test his theory, London calculated the transition 
temperature and heat capacity of the BEC (left figure). He 
showed that there is some agreement with liquid helium-4 
(right figure), at least in trend. 
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Left: Enss and Hunklinger, Low Temperature Physics, page 8, 
2005. 

Right: London, Physical Review, volume 54, page 947, 1938 
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BEC heat capacity. 



We shall look at a simplified version of London's treatment to 
calculate the BEC heat capacity. 

This involves finding: 

1. the temperature at which condensation starts, 

2. the change in number of atoms with temperature, and 

3. the heat capacity below condensation temperature. 
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The density of states 



We have seen a number of different density of states. 

For electrons, we need to multiply by 2 because of the spin 
states. For phonons, we multiply by 3 for the three possible 
polarisations. For photons, we multiply by 2 for the two 
polarisation states. 

So what do we do for atoms that are bosons? 

The answer is: we use the same density of state for the ideal 
gas: 

9(e) = — p— (2ms) 1 / 2 

This is the very first one that we have seen, before we have to 
include the additional effects of spin and polarisation. 
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When do the bosons condense? 



Using the Bose-Einstein distribution that we have learnt, we can 
calculate the temperature at which bosons condense into BEC. 

Recall the expression we first derived for the Bose-Einstein 
distribution: 

Hi 1 



Qi exp(-Ai - \ 2 Si) - 1 
When applying this to phonons, Ai was left out because the 
number of phonons is not fixed. 

Now that we are dealing with atoms, which are real bosons, the 
number is fixed and we have to keep Ai. This is usually written 
in terms of \i as follows: 

g(e)ds 



n(s)de = 



exp((£-/i)//c s T)-l 
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The chemical potential 



li is called the chemical potential. It is determined by the total 
number of particles. We get this by integrating: 

Jo exp(0 - ii)/k B T) - 1 

I have assumed that the ground state energy is zero, which is 
why I integrated from zero. 

This is not easy to solve. Instead, we shall make use of some 
approximations at low temperature. 

For a start, note that jj, could depend on temperature, since 
the above equation contains T. 
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The occupation number 



We focus on the occupation number in the integral: 



exp((e-fjL)/k B T) - 1 
First, note that \i cannot be positive, or else for small energy 
the exponential function would be less than one. Then the 
denominator would be negative. This means negative 
occupation number, which is not physical. 

Next, we know that when temperature is low enough, a large 
number of atoms would go into the ground state e = 0. The 
number would be as large as N, which is of the order of 1 mole 
(« 10 24 ). 

This means that /(0) would be very large, so that \i must be 
very close to zero. Then the exponential function would be 
close to 1 and the denominator would be close to zero. 
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The condensation temperature 

The formula for N ex is only able to compute the number of 
particles above the ground state. This is because the density of 
states g(e) in the integral is zero when s = 0. 



When a substantial fraction of the particles start going into the 
ground state, this has to be added separately. So N ex is the 
number of atoms in the excited states. Hence the subscript ex. 



We know that as soon as N ex becomes less than the original 
number N, then (N — N ex ) atoms start going into the ground 
state. 



Therefore, condensation takes place when N ex = N. 
Substituting this and solving for T, we get 

h 2 ( N \ 2 / 3 



T BE — 



2n;mk B V2.612V. 
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Number of atoms in the excited state 



So at very low temperature, it is safe to assume that \± = 0. 
We can then write 

ex Jo exp(s/k B T) - 1 

I have added a subscript ex to N. The reason will become clear 
in a moment. First, integrate this with the help of the Table of 
Integrals: 

_ f2-Kmk B T\ 3 / 2 



N, 



ex 



2.612V 

h z J 

Notice that T is in the denominator. This suggests that the 
total number of atoms decreases with temperature - that they 
are disappearing ! 



In fact, the above integral includes only atoms above the 
ground state e = 0. The "missing" atoms are in fact going into 
the ground state. 
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Transition temperature 



Substituting the mass of helium atom and the molar volume 
into the formula for the condensation temperature T BEl we find 
3.13 K. 

This is fairly close to the lambda point temperature of 2.18 K, 
at which liquid helium-4 changes to a superfluid, 
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and provides some support for the idea that the superfluid is a 
BEC. 
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Number of atoms in the ground state 



Above T BEl there are very little atoms in the ground state. 

Below T BE , this number increases until all atoms fall into the 
ground state. 



N 






^^^niut/tJ 2 ] 








1 



0 



Glazer and Wark, Statistical mechanics, p. 101 
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Below condensation temperature. 

The heat capacity below the condensation temperature is then 
obtained by differentiating with respect to T: 




Specific heat of an ideal Bose- Einstein gas. 



Note that the peak value is 1.926k B N. This is obtained by 
setting T = T BEl when N ex = N (all atoms are excited). 
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Heat capacity. 



Next, to calculate the heat capacity, recall the expression for 
the number of excited particles: 

g(s)ds 



'OO 



f 0 exp(s/k B T) - 1 

To find the heat capacity, we first find the energy U . The above 
integral is a sum over particle number in every energy interval. 



To find U, we need to multiply by the energy at each interval: 

eg(s)de 



'OO 



u = 



f 0 exp(e/k B T) - 1 
Integrating using the table of integrals gives 

T 5/2 



U = 0.7704k B N 



T 



3/2 



BE 
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Above condensation temperature. 



Above condensation temperature T BE , we cannot use the same 
formula for the heat capacity. This is because it is derived 
assuming that the chemical potential ji w 0, which is only true 
when condensation starts. 

Above T BEl ji changes. Albert Einstein has published a formula 
for this in 1924 (right half of curve): 




Specific heat of an ideal Bose-Einstein gas. 



At high temperatures, the heat capacity reaches the ideal gas 
value of 3Nk B /2. 
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The search for condensate. 



After London's work, the search for BEC in liquid helium 
started. Calculations and experiments were carried to 
determine if BEC really existed in the superfluid. 

Measuring the amount of BEC is difficult. The method is to 
use neutron scattering to determine the velocities of the helium 
atoms. 

In a BEC, the helium atoms are in the ground state - so their 
velocities are all zero. If the measurement can show that a 
significant fraction of the superfluid contains atoms with zero 
velocity, then we have found the condensate. 
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Helium-4 condensate 

The figure shows the theory and experimental results in the 
1970s and 1980s. The condensate has indeed be found. 
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Left: Tilley and Tilley, Superfluidity and Superconductivity, 
(2003) p. 64. 

However, even at 0 K, the condensate is only 10% of the 
superfluid. The reason is attributed to the interaction between 
atoms - it can only be 100% if the helium is an ideal gas at 0 K. 
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Neutron scattering. 



A quick look at how neutrons can be used to determine the 
velocities of atoms would be helpful. 

We cannot see the movement of the atoms in the superfluid, 
but we can measure the energy and direction of neutrons 
scattered from the superfluid. 













► 





Using momentum conservation, we can deduce the velocities of 
the helium atoms that the neutrons collided with. 
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Helium-4 superfluid 



Recall that London have proposed BEC as an explanation for 
superfluidity, because it behaves as a single wavefunction. 

If BEC is only 10% of the superfluid, that means only 10% of 
the superfluid is a single wavefunction. 

So why does 100% of liquid helium-4 behave as a superfluid at 
0 K? 

This problem is not fully understood until today. 
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Two-fluid model. 



Even without a full understanding, it is still possible to make 
models that are useful for predicting the behaviour. This would 
in turn provide some understanding. 

One month after London proposed the BEC theory for liquid 
helium in 1938, Laszlo Tisza proposed a two-fluid model for 
the superfluid. 




The idea is that below the 2.18 K transition temperature, 
normal and superfluid coexist. As temperature falls, the normal 
fluid decreases. At 0 K, the superfluid increases to 100%. 
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Andronikashvili's experiment. 



"According to this model, helium II behaves as if it were a 
mixture of two completely interpenetrating fluids with different 
properties, although in reality this is not the case. 

To avoid any misunderstanding, it must be clearly stated at the 
outset that the two fluids cannot be physically separated; it is 
not permissible even to regard some atoms as belonging to the 
normal fluid and the remainder to the superfluid component, 
since all 4He atoms are identical." 
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Connection with BEC. 



The two-fluid model is closely related to the BEC idea. Recall 
the left figure, showing the condensate going to 100% as 
temperature falls from T BE to 0 K. 




Right: http : //www . yutopian . com/Yuan/TFM . html 

If we include the graph for the excited (not condensate) 
fraction, it would look like the figure on the right. This is just 
what the two-fluid model would suggest. 
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Two-fluid model. 

In this model, the fluids are supposed to have the following 
properties: 





Normal fluid 


Superfluid 


Density 




P$ 


Velocity 


V N 


v s 


Viscosity 


Vn 


Vs = 0 


Entropy 




S s = 0 



Guenault, Basic Superfluids (2003) p. 29 

Note that we are not supposed to think that we can physically 
separate normal fluid from superfluid, since the helium-4 atoms 
are not distinguishable. 

This makes it look rather unrealistic. However, following these 
assumptions, it has been possible both to explain many 
experiments. 
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Andronikashvili's experiment. 



In 1948, Elepter Andronikashvili measured the fraction of 
normal fluid in liquid helium using a stack of rotating discs 
(left). The small spacings between discs means that any normal 
fluid would be dragged along and slow down the oscillation. 




1.5 2.0 2.5 3.0 3.5 

Andronikashvili's experiment. Temperature 77 k 



Right: Lifshitz and Andronikashvili, A Supplement to Helium 
(1959) 

In the experiment, the oscillation period increases a lot below 
2.2 K. By attributing this to a decrease in normal fluid in the 
spacing, the fraction of normal fluid can be deduced (right). 
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Boiling helium. 



When pressure is reduced below vapour pressure, superfluid 
helium boils (even if temperature is below 2.2 K). 

When a liquid like water boils, we expect to see bubbles 
forming the liquid moving. 

When superfluid helium boils, we cannot see anything 
happening! 

The helium atoms leave the surface of the liquid. Somehow, 
there is no bubble and no movement at all. 
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Viscosity measurement. 



Viscosity can be measured by allowing liquid helium to flow 
through a very narrow capillary tube. From the velocity of flow, 
the viscosity can be calculated. The result is shown in this 
graph: 




12 3 4 
Temperature 77 K 



J.F. Allen, A.D. Misener, Proc. R. Soc. A172, 467 (1939) 
K.R. Atkins, Philos. Mag. Supp. 1, 169 (1952) 



Notice that viscosity drops to zero below 2.2 K 
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Fountain effect. 



In the setup below (left), liquid helium in a flask and in a bath 
is separated by compressed powder which is porous. 





A 




c, 










a 

Heater 










^Powder 




Left: Enss and Hunklinger, Low-Temperature Physics (2005) 
p. 21. Right: J.F. Allen, photographed 1971, unpublished 



The temperature is about 1 K. When the liquid in the flask is 
warmed very slightly with a heater, it shoots out of the flask 
like a fountain (right). 
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Viscosity experiments explained. 



We have seen 2 viscosity experiments. 

In the rotating discs experiment by Andronikashvili, viscosity 
decreases gradually as temperature falls below 2.2 K. This 
shows that the superfluid fraction increases gradually. 

In the capillary tube experiment, the viscosity drops suddenly to 
zero. This suggests that the liquid suddenly becomes 100% 
superfluid. 

The reason is that the capillary tube is too narrow. The normal 
fluid would experience too much viscosity trying to go through. 
So only the superfluid fraction passes through, and this 
experiences zero viscosity. 
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Silent boiling explained. 



When water boils, there is plenty of bubbles and movements. 
This happens because certain spots in the water gets a lot 
hotter than other parts. Boiling then starts a these hotspots, 
where bubbles are formed. 

If bubbles are not formed in boiling helium, it must be because 
there are no hotspots. This can only be possible if the liquid 
has extremely good thermal conductivity. Any heat around 
potential hotspots must be conducted away immediately. 
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Fountain effect explained. 



In the fountain experiment, when the liquid helium in the flask 
is warmed, the liquid shoots out like a fountain. 

Using the two-fluid model, when it is warmed, the superfluid 
fraction changes to normal fluid. There is now a difference in 
concentration of superfluid between the flask and the helium 
bath below. 

Since they are separated by compacted powder, superfluid can 
pass. Normal fluid cannot because of viscosity. 

Since the concentration of superfluid in the flask decreased, 
superfluid will diffuse through the pores in the compacted 
powder from the bath. The volume of liquid in the fask then 
increases quickly and shoots out like a fountain. 
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Silent boiling explained. 



If any spot is hotter than the surrounding, it could convert the 
nearby superfluid into normal fluid. The concentration of 
superfluid decreases. Superfluid from further away flows in to 
the hotspot, and the normal fluid flows away. 

This process is very efficient because the superfluid has no 
viscosity. The flow is limited only by the viscosity of the normal 
fraction. 

As a result, heat conduction is extremely fast, and no hotspot 
can form. This is why no bubble is produced and the superfluid 
helium boils silently. 
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Landau critical velocity. 



Landau's dispersion curve is used to deduce a critical velocity 
below which liquid helium is superfluid. 




'0 0,2 0.4 06 0.8 1,0 L2 1.4 1.6 1.8 ZO ZZ Zfl 2.6 2B 3.0 
MOMENTUM CHANGE 0 (A" 1 ) 



Henshaw and Woods, Physical Review, volume 121, p. 1266 (1961) 

The above graph is measured using neutron scattering. Let us 
now see if the critical velocity calculated from this is indeed 
real. 
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Critical velocity. 



The drag measured is plotted against velocity of the ion 




Allum, et al, Philosophical Transactions of the Royal Society A, vol. 284 (1977), p. 179 



The result shows that the critical velocity is about 45 m/s. 
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Experiment. 



The idea is to accelerate a body through superfluid He-4, and 
see if it experiences a drag when the critical velocity is reached. 




Allum, et al, Philosophical Transactions of the Royal Society A, vol. 284 
(1977), p. 179 

The idea is to inject a helium ion into the superfluid. This ion 
would attract surrounding atoms into a ball with about 100 
atoms. It is then accelerated using an electric field. 
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From the measured dispersion, 




MOMENTUM CHANGE 0 (A" 1 ) 

the following quantities are obtained: 

A 

excitation: — = 8.65 K 

k B 

Po 

momentum: — = 19.1 nm 

h 
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The Landau critical velocity is then given by 

A _i 
v L = — = 58 m s . 

Po 

This is fairly close to the velocity of 45 m/s measured using 
the helium ion, and provides support for Landau's theory. 
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Phase of a wavefunction 



In this simple wavefunction, 

^ = e ikx , 

kx is the phase. In other wavefunctions, the phase can depend 
in a more complex way on the position x. Let the phase be 
00). 

The phase is useful to help explain the behaviour of the 
superfluid. 

First, note that the simple wavefunction, as we move along x, 
the phase changes. A change of 2tt corresponds to one 
wavelength. 

Next, let us redo the calculation on the previous slide in terms 
in terms of the phase. 
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Rotational motion. 



We now look at what the idea of a macroscopic wavefunction 
means for rotational motion in a superfluid. 

Consider a simple case of uniform motion in one direction. The 
wavefunction may be represented by a plane wave: 

i/> = e ikx . 

To see that the velocity is indeed constant, we use the 
equation in quantum mechanics for finding momentum: 

ax 

Substituting the wavefunction and differentiating, we find 

p x = hk. 

This is the familiar relation for momentum in quantum 
mechanics. Dividing by mass of the helium atom m gives the 
velocity: 

hk 



v = 



m 
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First, replace the phase in e lkx by the more general function 

cj)(x): 

Then, applying the equation in quantum mechanics for 
momentum: 

—in — = pip. 
dx 

we find 

dcf> 

p = a — . 

dx 

Dividing by mass of the helium atom m gives the velocity: 

H d<p 
v = . 

m dx 

Next, consider what happens to the phase in a circular motion 



Liquid Helium-4 



39 



Rotational motion. 



Suppose the superfluid is flowing round in a circle. Imagine 
following a path along this circle. When we return to the same 
point, the wavefuncton 

must return to the same value. Then the phase 0 must either 
change by 0, or by a multiple of 2tv. 

Let us now apply this idea to the velocity equation: 

H dcj) 

v = . 

mdx 

When we integrate along the circular path, we get 

f v.dl = — A0 

J m 

where A0 is the change in phase, and dl is like a very short 
vector along the path. 
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Rotation is quantised. 

In this equation 

/ v.dl = n— 

J m 

the left hand side is called the circulation. Since n is an integer, 
this means that the circulation is quantised. 

Note the circulation is a measure of the amount of rotation. To 
see this, start with the simplest case when there is no rotation. 

Then the velocity v is uniform, with the same magnitude and 
direction in the whole fluid. Integrating along a circle, we find 
that half the time, we move in the direction of flow, the other 
half we move opposite to the flow. Since there is a dot product 
in the circulation integral, that means the integral would be 
zero. 

Next, consider uniform rotation of the superfluid. When we 
integrate along a circle, v and dl are always in the same 
direction. Then the integral would be nonzero. Also, if v is 
larger, then the integral is larger. 
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In this equation, 

h 



[ v.dl = — A0. 

J m 



m 

the dot product is used on the left hand side. This is because 
we can think of the wavefunction like waves on the sea. 

Only moving along direction v of the wave causes phase 
change. Moving along the phase front perpendicular to v gives 
no change in phase. The dot product picks out the component 
along v. 

We have seen that A</> must be 0 or a multiple of 2tt, so the 
equation becomes 

f 71 h ^ 

/ v.dl = — 2n7r = n — , 
J ra m 

where n is an integer. 
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Quantised vortex. 



When water rotates in a bucket, is looks like this: 




We call this a vortex. It is difficult to imagine a vortex like this 
can only rotate at certain speeds. Yet, this is exactly what 
London predicted with the wavefunction, if we have a bucket of 
superfluid helium. 

Subsequently experiments have shown that this is indeed true. 



Liquid Helium-4 



43 



Measuring vortices. 



One way to measure circulation is to look at the vibration of a 
wire in a superfluid. When the circulation changes, the 
vibration changes. Using this method, the circulation is plotted 
against angular velocity: 
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Karn, et al, Physical Review B, vol. 21 (1980), p. 1797. 

The distinct steps show that the circulation can only take 
discrete values, i.e. it is quantised. 
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Exercises 



Exercise 1 

Bose-Einstein condensation is achieved with a gas of rubidium 
atoms, which have an atomic weight of 85.47. The 
condensation temperature is 1.7 x 10 -7 K. What was the 
number density of the condensate? 

[ Atomic mass unit u is 1.6605 x 10 -27 kg.] 
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What we have learnt so far. 



1. Below 2.2 K, liquid helium-4 loses all viscosity. This new 
phase is called a superfluid. 

2. London used Bose-Einstein condensation to explain this. 
He calculated a transition temperature of 3.13 K. 

3. Theory and experimental search show that a condensate 
really does exist in superfluid helium-4, but only at 10%. 

4. The two-fluid model can be used to explain observations 
like silent boiling, fountain effect, and zero viscosity. 

5. Landau's prediction of a critical velocity is directly verified 
by experiments of moving ions. 

6. Rotational motion in superfluid helium is quantised. 
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Exercises 



Recall the formula for the condensation temperature: 

h 2 / N \ 3 / 2 

BE ~ 27rmk B V2.612W 
We can rearrange this to get the number density: 

V V h 2 J 

The condensation temperature is given as 

T BE = 1.7 x 10" 7 K. 

The mass is given by 

m = 85A7u = 85.47 x (1.6605 x 10" 27 ) kg. 

Substituting these into the above equation for the number 
density, we get 2.74 x 10 19 m -3 . 
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Exercises 



Exercise 2 

Show that below the condensation temperature, T BE , the heat 
capacity of a gas obeying Bose-Einstein statistics is given by 

\ 3/2 

C v = 1.93Nk B [-^- 

. BE. 



[You are given that 



roc x 3 l 2 

/ dx = 1.78. 

f 0 e x - 1 
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Exercises 



Define 



s 

X = 



k B T 

and substitute into the integral for energy: 

AmirV r- — . _ n r/o f°° x 3 / 2 dx 

We can now use the given result: 

roc x 3 ! 2 

dx = 1.78. 



Jo e x — 1 
This gives: 

u = l!^^ (fcBT )5/2 x 17Q 

The heat capacity can then be obtained by differentating this 
with respect to temperature. First, we express this in terms of 
the condensation temperature. 
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Exercises 



Recall that at very low temperatures, the chemical potential \i 
is approximately zero, and the number of particles above the 
ground state is given by: 

roc g(s)ds 



N ex = 



'o exp(s/k B T) - 1 
To obtain the energy, we insert e into the integral 

Jo exp(s/k B T) - 1' 
Substituting the density of states 

AmnV 



gives 



g{e) = — p— (2ms) 1 / 2 



_ 4-mitV f°o ^J2rne 3 / 2 de 



h 3 Jo exp(e/k B T) - 1 
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Exercises 



The condensation temperature is given by: 

_ h 2 / N \ 3 / 2 
TBE ~ 2irmk B VZ612V J 
Make V the subject 

v = — ( ) 3/2 

2.61 \2ixmk B T BE ) 
and substitute into the above result for U: 

2 1.78 T 5 / 2 
U = -= Nk B —^ 

Then differentiating with respect to T: 

T \3/2 



C V = 1.93Nk B 



■ 



T BE, 

we get the heat capacity for the Bose-Einstein condensate. 
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